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Semiclassical coupled-channels equations. Consider H = Hq + V, where H is the 
Hamiltonian composed by a non-perturbed Hamiltonian Hq and a small perturbation 
V. The Hamiltonian Hq satisfies an eigenvalue equation Hq^ u = E n y m whose eigen- 
functions form a complete basis in which the total wavefunction that obeys the 
Schrodinger equation H*¥ = ihdm Jdt, can be expanded: 

V = Y,an(t)Vne- iEnt/h . (1) 

Inserting this expansion in the Schrodinger equation, one obtains 

ih^dnVne- 1 ^'* = Yya n ^ n e- lE "^\ (2) 

n n 

with d n = da n (t)/dt. Using the orthogonalization properties of the we multiply §2§ 
by y/f and integrate over the coordinate space to get the coupled-channels equations 

<*k(t) = -jrEa«(0 Vkn(t) f^ 1 , (3) 

n n 

where the matrix element (dx is the coordinate volume element) is Vk n = j ^Vy/ n dx. 

Often, the perturbation V is very small and the system, initially in state n = 0, does not 
change appreciably. Thus we can insert a n = 8 n o in the right-hand side of the equation 
(|3]), which yields 

a k (t)= f dtV kn (t) ^*-W*. (4) 

This is the first-order perturbation theory result. Obtaining the coefficients a n , the 
probability of occupation of state n is simply given by \a n {°°) | 2 . Isn't that fun? So simple, 
yet so powerful! In energy space, things are not so different (I will be back to this later). 

Multipole expansions. Let us now consider a point particle with charge Z p e at a 
distance r from the center of a charge distribution where we put the origin of our frame 



of reference. At a point r' from the origin let the charge density be p(r'). The potential 
created by the charge Z p e located at r, averaged over the charge distribution, is 

cW p i M ~ + 7 + l^ + '''' () 

where p = / r'p(r')d 3 r' and = / (3r-r^- — r' 2 )p(r')d 3 r' are the dipole and quadrupole 
moments of the charge distribution, respectively. In the last step we have expanded the 
factor l/|r — r'| for r <C /, rearranged terms, and used / p(r')d 3 r' = Zje. No big deal. 
Try deriving it yourself. I bet you can do it even after eating pasta with chianti. 

Semiclassical what?!. In the semiclassical approximation, one assumes that the 
projectile's coordinate r can be replaced by r(t), following the classical trajectory of the 
projectile. The dependence on r' is used to treat p and Qu as operators. The excitation of 
nucleus T by nucleus P is obtained by using the matrix element (/|Vc(r',f)|i), where |i) 
( |/)) is initial (final) state of nucleus T. In layperson's terms, it means the replacement 
of the static density p(r') by the transition density Pfi(r r ) = ^Py-(r / )^/(r / ) in Eq. (|5]), 
where Q¥f) is the initial (final) wavefunction of nucleus T. 

The validity of the semiclassical approximation relies on the smallness of the wave- 
length jk = X/2n of the projectile's motion as compared to the distance of closest ap- 
proach between the nuclei. Let us consider a central collision. The distance of closest 
approach is a = 2ao = ZpZje 1 /E, where E = mv 2 /2 is the kinetic energy of relative 
motion between the projectile and the target and m is the reduced mass. We can use 
the so-called Sommerfeld parameter, r\ = a$j /i to measure the validity of the semi- 
classical approximation. Since =fi/p = h/mv, we have r\ = ZjZpe 1 / mv ^> 1 as the 
condition for validity of the semiclassical approximation. It is easy to verify that for 
nucleus-nucleus collisions this condition is valid for most cases of interest. In summary, 
the "semi" from semiclassical means "quantum". The "classical" means that the scatter- 
ing part is treated in classical terms. The later is well justified in most situations. But 
don't worry, we can also do all quantum easily, as I show later. 

Low energy central collisions. The fun part starts here. Consider the transition of 
the ground state J = of a deformed nucleus to an excited state with J = 2 as a result of a 
frontal collision with scattering angle of = 180°. From Eq. ([5]) the perturbing potential 
is V = jZ p e 2 Qif /r 3 . According to Eq. Q, the excitation amplitude to first order is 

Z P e 2 Q if fe im 

At an scattering of 9 = 180° a relationship exists between the separation r, the velocity 

1/2 

v, the initial velocity vo, given by v = drj dt = ±vo (l—a/r) 1 (show this using energy 
conservation). If the excitation energy is small, we can assume that the factor e mt in (6j) 
does not vary much during the time that the projectile is close to the nucleus. The orbital 
integral is then solved easily by the substitution u = l—a/r, resulting in (you can do 
this integral, I know) 

_ 4Z p e 2 Q if _ AQjfE 2 
3ihvoa 2 3Z p e 2 fiVQZj 



The differential cross section is given by the product of the Rutherford differential cross 
section at 180° and the excitation probability along the trajectory, measured by the 
square of a,y. That is, do/da\ e=mo = dOR/dQ. |e=i80° x \ a if\ 2 - One obtains 



do 
dQ. 



m E\Qi f \ 
0=180° lSH 2 Zi 



(8) 



where mo is the reduced mass of the projectile+target system. Oops! This expression is 
independent of the charge of the projectile. But you have always heard that heavy ions 
(large Zp) are more efficient for Coulomb excitation. What is wrong? Show that there is 
nothing wrong. This formula is right and what you've heard is correct. 

General multipole expansion. Instead of Eq. (|5]), for electric excitations an exact 
multipole expansion can be carried out with the help of spherical harmonics: 

A-rr 1 

V(r) = Z-—^ TI Yl M (r)^(EL,M). (9) 

LM 

The electric multipole moment of rank L = 0, 1 , • • • is given by ^(EL,M) = 
fd 3 r'p(r'y L Y LM (r'), where M = -L,-L+l,.. .,+L. 

In semiclassical calculations, r is replaced by a time-dependent coordinate along a 
Rutherford trajectory for the relative motion, r(t) . As before, we can proceed to calculate 
the excitation amplitude for the transition of the target from a state with energy E{ to a 
state with energy Ef by replacing p by p,y (or M by ^/) and integrating over time 
from — oo to +oo. The cross section is obtained by squaring the excitation amplitude, 
summing it over final and averaging over the initial intrinsic angular momenta of the 
target. Oh, don't forget to multiply it by the Rutherford trajectory. One gets (with some 
omitted factors) 



da 



Z Z P B(EL) iw<»/ ; )r, (10) 



where CO/,- = (Ef — Ei)/h and B(EL) ~ / dr' r' L pi(r') is the reduced transition prob- 
ability. By the way, Pl(^) is the radial part of the L-pole component of the transition 
density p i /(r / ). The factor I E l involves a sum over M of the orbital integrals lELM(tt>fi) 
given by 

hLM{(Ofi) = jdt - / ^YufLM ir(t))e i(0fit . (11) 



Virtual photon numbers. Integration of ( 10 ) over all energy transfers Ey = Ef — Ej, 
yields 

d^rd^^rdEydnEE {E)a EL {E) (12) 

da da {tij Ej da K v 7 K yh 



where Oy L are the photonuclear cross sections, 



Oy L (Ey) ~ Ey L ~ l B(EL) 



(13) 
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FIGURE 1. Left: Equivalent photon numbers per unit area incident on 208 Pb, in a collision with 16 at 
100 MeV/nucleon and with impact parameter b = 15 fm, as a function of the photon energy E = hco. The 
curves for the El, E2 and Ml multipolarities are shown. Right: Total number of virtual photons for the 
El multipolarity, "as seen" by a projectile passing by a lead target at impact parameters b m i„ = 12.3 fm 
and larger (i.e., integrated over impact parameters), for three typical bombarding energies. 



The virtual photon numbers, nEL(Ey), are given by 

(E 7 )~zl\I EL (E r )\ 2 . (14) 



dn E L /T7 , 7 2 



The dependence of the cross section on the deflection angle is included in dnEL/dO.. 
Since for a Rutherford trajectory the deflection angle is related to the impact parameter 
by b = ciq cot 6/2, we can write the cross section in terms of an impact parameter 
dependence by using dn£L/27lbdb ~ sm 4 (6)dnEL/dCl. 

The formalism above can be extended to treat magnetic multipole transitions, ML. It 
is much more complicated (involves currents, spins): hard work, but straight-forward. 

The reactions induced by real photons include the contribution of all multipolarities 

with the same weight, i.e., Oy{Ey) = L\E/M.L G y^ M,L {Ey), whereas according to Eq. ( fT2] ), 



the excitation by virtual photons (e.g. Coulomb excitation) has different weights, n E j M , 
for different multipolarities. 

Figure ([T]) (left) shows the equivalent photon numbers per unit area, dnEL/^bdb, 
incident on^ 08 Pb, in a collision with 16 at 100 MeV/nucleon and with impact parameter 
b=\5 fm, as a function of the photon energy E = h(0. The curves for the El, E2 and Ml 
multipolarities are shown. One sees that there is a cutoff for excitation energies beyond 
the adiabatic limit , i.e., Ey < yfiv/b. On the right we show the total number of virtual 
photons, n E L — ft? . dbdnEhl (2nbdb) , for the E 1 multipolarity, "as seen" by a projectile 
passing by a lead target at impact parameters b m i„ = 12.3 fm and larger (i.e., integrated 
over impact parameters), for three typical bombarding energies. Lesson: more and more 
high energy photons are available as the beam energy increases. 

We can easily understand the origin of the adiabatic condition by investigating the 



orbital integral, Eq. ( fTTj ). Notice that for times larger than t exc =1/(0 the integral 



oscillates too fast and Ielm is small. For collisions at low energies, the collision time 
is given by t co u = ciq/v, where ciq = ZpZje 2 /2E CM _. Thus, the excitation is possible if 
tcoll/texc < 1 otherwise the system will respond adiabatically (i.e. nothing interesting 
happens). This condition is called the adiabatic condition, coao/v < 1. For collisions 
at high energies, nuclei follow nearly straight-line orbits and it is more appropriate 
to use the impact parameter, b, as a measure of the distance of closest approach. The 
collision time is t co u = b/yv, where 7=1/ yl — v 2 /c 2 is the Lorentz contraction factor. 
Thus the adiabatic condition becomes cob/yv < 1. Assuming an impact parameter of 
20 fm, states with energy up to yhv/(20 fm) can be appreciably excited. Thus, even 
for moderate values of y, i.e., y = 1 — 2, it is possible to excite giant resonances. With 
increasing bombarding energies, ultraperipheral collisions can access the quasi-deuteron 
effect, produce deltas, mesons (e.g., JW), even the Higgs boson. Whatever! You name 
it. 



Nuclear response to multipolarities. The response function is 

2 



B(EL) 



Jr L p if d 3 r 



(15) 



where p,y = *PpP,- is the transition density. A simple estimate can be done for the 
excitation of high multipolarities by assuming that the wavef unctions have the form *Pj = 
which yields B(EL) ~ R 2L , or, from Eq. ([13]), a[ ~ (kR) 2L , where k = 



Ey/hc. Thus, (Tl+i/ctl ~ (kR) . Usually kR <C 1 (long-wavelength approximation) for 
low-lying states and we see that the cross sections decrease strongly with multipolarity. 

It is useful to estimate the total photoabsorption cross section summed over all tran- 
sitions \i) — ► \f). Such estimates are given by sum rules (SR) which approximately 



ZfiEf-Et) 



i) 



Here the 



determine quantities of the following type: 

transition probabilities for an arbitrary operator are weighted with the transition en- 
ergy. For such energy-weighted sum rules (EWSR), , a reasonable estimate can be 
derived for many operators under certain assumptions about the interactions in the sys- 
tem. Using the completeness of the intermediate states, the commutation relations be- 
tween the Hamiltonian and the operator &, assuming that the Hamiltonian does not 
contain momentum-dependent interactions, one gets the EWSR for the dipole operator, 
d z = rFio(f), y[d z ] = ^ a ^ 2 e 2 /2m a , where the sum extends over all particles with mass 
m a and charge e a .This is the old Thomas-Reiche-Kuhn (TRK) dipole SR. 

We have to exclude the center-of-mass motion. Therefore our z-coordinates should 
be intrinsic coordinates, z a =>• z a —R z , where R z = L fl z a /A. Hence, the intrinsic dipole 

moment is d z = £ a e a (z a -R z )=e £ p z p - (Ze/A) (^£ p z p + E„ z n J ■ This operator can be 

rewritten as d z = e„ £ p z p + e n z n where protons and neutrons carry effective charges 



eN/A, 



-eZ/A. (Weird, no? Think about it.) This yields the dipole EWSR 



y[d z ]=Y,E fi \dz fi \ 
f 



h 2 e 2 NZ 
2myy A ' 



(16) 
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FIGURE 2. Left: Coulomb excitation cross section of giant resonances in 40 Ca projectiles hitting a 238 U 
target as a function of the laboratory energy per nucleon. The dashed line corresponds to the excitation of 
the giant electric dipole resonance, the dotted to the electric quadrupole, and the lower line to the magnetic 
dipole. The solid curve is the sum of these contributions. Right: Virtual photon numbers for the electric 
dipole multipolarity generated by 84A MeV 17 O projectiles incident on 208 Pb, as a function of the center- 
of-mass scattering angle. The solid curve is a semiclassical calculation. The dashed and dotted curves are 
eikonal calculations with and without relativistic corrections, respectively. 



where is the nucleon mass. The factor (NZ/A) is connected to the reduced mass for 
relative motion of neutrons against protons as required at the fixed center of mass. 



The EWSR ( |16| ) is what we need to evaluate the sum of dipole cross sections for real 
photons over all possible final states |/). Taking the photon polarization vector along the 
z-axis, we obtain the total dipole photoabsorption cross section 

oL^fdE^lTC 2 ^^. (17) 



f 



This universal prediction on average agrees well with experiments in spite of crude- 
ness of approximations made in the derivation. One should remember that it includes 
only dipole absorption. For the E2 isoscalar giant quadrupole resonances one can derive 
the approximate sum rule / dE r ojjQ R (E 7 ) /E* ~ 0.22ZA 2 / 3 jih MeV -1 . 

Resonances. A simple estimate of Coulomb excitation of giant resonances based 
on sum rules can be made by assuming that the virtual photon numbers vary slowly 
compared to the photonuclear cross sections around the resonance peak. Then 

° C ~ ^EgdR^ f dEr °GDR( E y) + n E2(E GQR )E GQR J ^^Q^Ey). (18) 

In figure we show the Coulomb excitation cross section of giant resonances in 40 Ca 
projectiles hitting a 238 U target as a function of the laboratory energy per nucleon. The 
dashed line corresponds to the excitation of the giant electric dipole resonance, the dotted 
to the electric quadrupole, and the lower line to the magnetic dipole which was also 



obtained using a sum-rule for Ml excitations 0. The solid curve is the sum of these 
contributions. The cross sections increase very rapidly to large values, which are already 
attained at intermediate energies (~ 100 MeV/nucleon). 

As with giant dipole resonances (GDR) in stable nuclei, one believes that pygmy 
resonances at energies close to the threshold are present in halo, or neutron-rich, nuclei. 
The hydrodynamical model predicts (31 for the width of the collective mode T = hv/R, 
where v is the average velocity of the nucleons inside the nucleus. This relation can be 
derived by assuming that the collective vibration is damped by the incoherent collisions 
of the nucleons with the walls of the nuclear potential well during the vibration cycles 
(piston model). Using v = 3\p/4, where vp = \/2Ep jm^ is the Fermi velocity, with 
Ef = 35 MeV and R = 6 fm, one gets T ~ 6 MeV. This is the typical energy width a 
giant dipole resonance state in a heavy nucleus. In the case of neutron-rich light nuclei 
v is not well defined. There are two average velocities: one for the nucleons in the core, 
v c , and another for the nucleons in the skin, or halo, of the nucleus, v^. Following Ref. 
[BM93], the width of momentum distributions of core fragments in knockout reactions, 
o t , is related to the Fermi velocity of halo nucleons by vf = a/ 5o 2 /m^. Using this 
expression with o c ~ 20 MeV/c, we get T ~ 1 MeV, in accordance with experiments. 
Usually such modes are studied with the random phase approximation (RPA). 

Eikonal waves. The free-particle wavefunction y ~ e lk r becomes "distorted" in 
the presence of a potential V(r). The distorted wave can be calculated numerically by 
performing a partial wave-expansion solving the Schrodinger equation for each partial 
wave, i.e., if y/ = Lim(Xi(r)/r)Y Im (r), then 



where 



EL 

dr 2 
h 2 



Xl{r) = 



E-V{r)- 



l{l+l)h 2 
2/lr 2 



1/2 



(19) 



(20) 



with the condition that asymptotically \j/(r ) behaves as a plane wave. 

The solution of ( [T9] ) involves a great numerical effort at large bombarding energies E . 
Fortunately, at large energies E a very useful approximation is valid when the excitation 
energies AE are much smaller than E and the nuclei (or nucleons) move in forward 
directions, i.e., <^ I . Calling r = (z, b), where z is the coordinate along the beam 
direction, we can assume that yA(r) = e ife 0(z,b), where is a slowly varying function 
of z and b , so that | V 2 | <^ k \ V0 1 . In cylindrical coordinates the Schrodinger equation 
for i// becomes 



2ike lk ^ + e^ 
oz 



dz 2 



v 2 h <t> 



2m 

J 2 



Ve ikz (j)=0 



or, neglecting the 2nd and 3rd terms, we get d(j)/dz= —iV(r)/hv 0, whose solution is 



(j) = exp 



i 



V(b,z')dz' 



(21) 



That is, 



yr(r) =exp{ikz + iz(p,z)}. 



(22) 



This is the eikonal function, where 

X(b,z) = -^-J^V(b,z')dz' (23) 

is the eikonal phase. Given V(r) one needs a single integral to determine the scattering 
wave. Do you have any idea how many people made their lives from the eikonal waves?. 
Well, don't ask, don't tell. By the way, some people call anything carrying an eikonal 
wavefucntion by "Glauber" theory. 

Quantum scattering. Defining r as the separation between the center of mass of the 
two nuclei and r' as the intrinsic coordinate of the target nucleus, the inelastic scattering 
amplitude to first-order is given by [0 

f(0) = ^r Jd 3 r rfV (< } (r) W) | *(r, r') | < +) (r) ( -(r')) , (24) 

where 4^, ^ (r) and <i>^ (r) are the incoming and outgoing distorted waves, respectively, 
and 0(rj is the intrinsic nuclear wavefunction of the target nucleus. Looks complicated. 
But that is the way we calculate quantum scattering amplitudes. Sometimes one calls 
this the Distorted Wave Born approximation (DWBA). 

At intermediate energies, AE/E[ ab <C 1, and forward angles, 6 < 1, we can use 
eikonal wavefunctions for the distorted waves. Corrections due to the extended nuclear 
charges can also be easily incorporated 10. The results can also be cast in the form of 



Eq. ( |T2[ ). Trust me on this one. 

In figure [2] we show the virtual photon numbers for the electric dipole multipolarity 
generated by 84A MeV 17 projectiles incident on 208 Pb, as a function of the center- 
of-mass scattering angle. The solid curve is a semiclassical calculation. The dashed and 
dotted curves are eikonal calculations with and without relativistic corrections, respec- 
tively (relativity? Well, no space to explain it here. Next school.). One sees that the 
diffraction effects arising from the quantum treatment of the scattering change consid- 
erably the differential cross sections. The corrections of relativity are also important. 
However, for small excitations both semiclassical and quantum scattering yield similar 
results for the differential cross section, as shown in figure [3] (left). 

Single particle and collective response. Assume a loosely-bound particle described 
by an Yukawa of the form exp(— r\r)/r, where r\ = ^2jlb c S/h, \i\, c is the reduced mass 
of (particle b + core c), S is the separation energy. This is a reasonable assumption for 
the deuteron and also for other neutron halo systems. We further assume that the final 
state is a plane-wave state (i.e., we neglect final state interactions) y/y = (q|r) = e' q r . 



The response functions for electric multipole transitions, calculated from Eq. ( |T5| ) is 

dB(EL;E 7 ) Vs(E y - S) L+1 / 2 



dE y E^ L + 2 



(25) 
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FIGURE 3. Left: Comparison of quantum (solid) and semiclassical (dashed) calculations of Coulomb 
excitation by E x =1.5 MeV of 8 B projectiles incident on lead at 50 and 250 MeV/nucleon, respectively. 
Right: Square of the integrals over coordinate, |=^J 2 , used in the calculation of the response function (T5\f 
in the single-particle model for a light, weakly-bound nucleus. Different scattering lengths and effective 
ranges were assumed for the final state. 



The maximum of this function occurs at Eq L = (L+ |)S/(L+ 1) ~ S. 

In figure [3] (right) we show a similar calculation as described above, but accounting for 
final state interactions in the form of scattering lengths and effective ranges (I wish I had 
more space to explain that, too.). The integrals over coordinate, denoted by y sp show a 
strong dependence on the final state interactions. The strong dependence of the response 
function on the effective range expansion parameters makes it an ideal tool to study the 
scattering properties of light nuclei which are of interest for nuclear astrophysics. 

The Coulomb dissociation method. As discussed above, the Coulomb breakup 
cross section for a + A — > b + c + A can be written as 



d -^^=N^(E Y ;e)o^ a ^ b+c (E r ), (26) 



where Ey is the energy transferred from the relative motion to the breakup, and 
°y+a -> b+c^Py) is the photo-dissociation cross section for the multipolarity TtL and 
photon energy Ey. Time reversal allows one to deduce the radiative capture cross section 
b + c -> a + 7 from a^ a ^ b+c (E Y ) , i.e., 



2(2./*+ 1) k 



Ob+c^ y+a ~ (2^ + l)(2j c +l) B° na "> b+ " (27) 

where ky = Ey/hc is the photon wavenumber, and k = (Ey — B)/h is the wavenum- 
ber for the relative motion of b+c. Except for the extreme case very close to the threshold 
(k — > 0), we have ky k, so that the phase space favors the photodisintegration cross 



section as compared to the radiative capture. Direct measurements of the photodisin- 
tegration near the break-up threshold do hardly provide experimental advantages and 
seem presently impracticable. On the other hand the copious source of virtual photons 
acting on a fast charged nuclear projectile when passing the Coulomb field of a (large Z) 
nucleus offers a way to study cross sections close to the breakup threshold. 

This method was introduced in Ref. [6] and has been tested successfully in a num- 
ber of reactions of interest to astrophysics. The most celebrated case is the reaction 
7 Be(/?, y) 8 B (see figure Q left). This reaction is important because it produces 8 B in the 
core of our sun. These nuclei decay by emitting high energy neutrinos which are one of 
the best probes of the sun's interior. The measurement of such neutrinos is very useful 
to test our theoretical solar models. 



Semiclassical CDCC. The Coulomb dissociation method is specially useful if first- 
order perturbation theory is valid. If not, one can still extract the electromagnetic matrix 
elements involved in radiative capture reactions. But a much more careful analysis of the 
high-order effects needs to be done. This can be accomplished by using a time-dependent 
discrete states are defined as 

| (j) b ) = e ~ iEb 'l h \b) , for bound states 

and =e~ iE J t/h J Tj(E) \E,JM) for continuum states, (28) 

where \E,JM > are continuum wavefunctions of the projectile fragments (with or with- 
out the interaction with the target), with good energy and angular momentum quantum 
numbers £, JM. The functions Tj(E) are assumed to be strongly peaked around an en- 
ergy Ej in the continuum. Therefore, the discrete character of the states \QjjAf) (together 
with \<j>b)) allows an easy implementation of the coupled- states calculations (see fig.|4j 



right). Calling them all together by \a), the orthogonality of the discrete states ( |28] ) is 
guaranteed if / dE T a (E) Tp(E) — 8 a p. Writing the time-dependent Schrodinger equa- 
tion for *F(f) = Y,a a a(t)<l>a, taking the scalar product with the basis states and using 
orthonormality relations, we get the coupled-channels equations ([3]). The problem of 
higher-order effects has been solved in this way for several cases. It is known as Semi- 
classical Continuum Discretized Coupled-Channels (S-CDCC) method. You can drop 
the "S" if you want. 

Schrodinger equation in a lattice. Another treatment of higher-order effects as- 
sumes solving the Schrodinger equation directly by discretizing space and time. This 
equation can be solved by a finite difference method assuming that the wavef unction 
can be expanded in several bound and unbound eigenstates \a), as before. A truncation 
on the sum is obviously needed. To simplify, we discuss the method for one-dimensional 
problems. The wave function m a at time t + At is obtained from the wave function at 
time t, according to the algorithm Q 



1 (0 \ At 



1 (0 \ At At 
IX 2hx hx 



(29) 
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FIGURE 4. Left: S-factor for the 1 Be(p,yfB reaction. The solid red circles are data obtained with the 
Coulomb dissociation method. The curves represent different theoretical methods for the reaction. Right: 
The reaction of a weakly-bound nucleus can lead to several steps of excitation in the continuum. 



In this equation z = HAt /4/i bx (Ax) 2 and SWcc^) = L (a \V\ a'} (t), with V being the 

a' 

time dependent potential, responsible for the transitions. Vq is part of Hq. 

The wave functions m a (x,t) are discretized in a mesh in space, with a mesh-size Ax. 

The second difference operator A( 2 ) is defined as = *pJ/ +1) (?) + ^iT^W - 

2^(0, with ^=V a (rj,t). 

The wave function calculated numerically at a very large time will not be influenced 
by the Coulomb field. The numerical integration can be stopped there. The continuum 
part of the wave function is extracted by means of the relation (and normalized to unity) 

V c (r,/) = \V-V gs < V gJ | >] [1- |< V gs | >\ 2 y l/2 (30) 

where is the initial wave function. This wave function can be projected onto an 
(intrinsic) continuum state to obtain the excitation probability of the state. 

Eikonal CDCC. To get quantum dynamical equations to treat higher-order effects, 
one discretizes the wavefunction in terms of the longitudinal center-of-mass momentum 
k z , using the ansatz 

< P = £^a(z,b) exp(iW) (3D 

a 

In this equation, (z,b) is the projectile's center-of-mass coordinate, with b equal to 
the impact parameter. (£) is the projectile intrinsic wavefunction and (k,K) is the 
projectile's center-of mass momentum with longitudinal momentum k and transverse 
momentum K. 



Neglecting terms of the form V 2 y a (z,b) relative to ikdz^a (z,b), the Schrodinger 
(or the Klein- Gordon) equation reduces to 

mv d^ a (z,b) = £( a \V\a') y d {z,b) A**-*")*. (32) 

° z a' 

These are the eikonal-CDCC equations (E-CDCC). They are much simpler to solve than 
the complicated low-energy CDCC equations because the z and b coordinates decouple 
and only the evolution on the z coordinate needs to be treated non-perturbatively. Of 
course, I lied and there are other complications (angular momentum coupling, etc.) 
hidden below the rug. If quantum field theorists can do it, why can't we? 

The matrix element (a \V\ <x') is Lorentz invariant. Boosting a volume element from 
the projectile to the laboratory frame means J 3 ^ — ► J 3 ^ /y. The intrinsic projectile wave- 
function is a scalar and transforms according to <j) a (% x , £ z ) — ► (j) a (% x , y^ z ), while 
V, treated as the time-like component of a four- vector, transforms as V (b,z; — > 
yV (b,z;t^x,^y,y% z ). Thus, redefining the integration variable z in the laboratory as 
£2 = y^z leads to the afore mentioned invariance. We can therefore calculate (oc\V\ a') 
in the projectile frame. Good Lord. That makes calculations so much easier. 

The longitudinal wavenumber hk a c ~ [E 2 —M 2 c 4 ) 1//2 also defines how much energy 
is gone into projectile excitation, since for small energy and momentum transfers k! a — 
k a ~ (E' a — E a ) jhv. In this limit, eq. ([32]) reduces to the semiclassical coupled-channels 
equations, Eq. Q, if one uses z = vt for a projectile moving along a straight-line classical 
trajectory, and changing to the notation y a (z,b) = a a (t,b), where a a {t,b) is the time- 
dependent excitation amplitude for a collision wit impact parameter b. Isn't that great! 
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